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A search for the Higgs Boson in H— t-ZZ— mode 
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A search for a Higgs boson in the decay channel H — > ZZ'*' with each Z boson decaying to an 
electron or muon pair is presented using pp collisions from the LHC at i/s — 7 TeV. The data 
analyzed correspond to an integrated luminosity of 1.13 ± 0.07 fb -1 recorded by the CMS detector 
in 2010 and 2011. The search covers Higgs boson mass (mn) hypotheses of 110 < mn < 600 GeV/c 2 . 
Fifteen events are observed, while 14.4 ± 0.6 events are expected from standard model background 
processes. Upper limits at 95% CL on the cross section x branching ratio for a Higgs boson with 
standard model-like decays exclude cross sections from about one to two times the expected standard 
model cross section for masses in the range 180 < mn < 420GeV/c 2 . Reinterpreted in the context 
of the standard model with four fermion families a Higgs boson with a mass in the range 138 -162 
GeV/c 2 and 178-502 GeV/c 2 is excluded at 95% CL. 

I. INTRODUCTION 

The standard model (SM) of electroweak interactions [T7] predicts the existence of a scalar boson, the Higgs 
boson, associated to the spontaneous electroweak symmetry breaking QJj]. The mass mn of this scalar boson is 
a free parameter of the theory. The inclusive production of SM Higgs bosons followed by the decay H— iTSZS*} is 
expected to be a main discovery channel at the CERN LHC pp collider. 

The direct searches for the SM Higgs boson at the LEP e + e~ collider have lead to a lower mass bound of 
tor > 114.4 GeV/c 2 (95% CL) [BJ. The direct searches by the DO and CDF experiments at the Tevatron exclude 
the mass range 158 < mjj < 173 GeV/c 2 (95% CL) [2J. Indirect constraints from precision measurements favour 
the mass range m H < 185 GeV/c 2 (95% CL) [JJ. 

A search for a SM Higgs boson in the four-lepton decay channel, H— >ZZ^*) — > l^Pl ± £ T with £,£' = e or /i, in 
short H— » A£, is presented. The analysis is designed for a Higgs boson mass in the range 110 < mjj < 600 GeV/c 2 
and uses data collected by the CMS experiment during 2010 and 2011 at the LHC collider, with pp collisions 
at y/s — 7 TeV. The sample of data corresponds to an integrated luminosity of £ = 1.13 ± 0.07 fb _1 . 

The sample of events with four reconstructed leptons contains an irreducible contribution from ZZ^ pro- 
duction via qq and gg fusion processes. Potential reducible background contributions are from Zbb and 
tt — > W + 6W _ 5, with the W undergoing a leptonic decay, where the final states contain two isolated lep- 
tons. Reconstructed 4£ events can also arise from instrumental background such as Z+jets or WZ + jet(s) or 
from the production of multiple jets in QCD hard interactions, where jets are misidcntified as leptons. 

A detailed description of the CMS detector can be found elsewhere [4]. 

The trigger used in this analysis evolved in response to the increasing instantaneous luminosity. The presence 
of a pair of electrons with Ef ri > 17 and Efp 2 > 8 GeV, or a pair of muons with j>t\ > ^ ano - Pt2 ■ > 
8GeV/c, is required for a major fraction of the data sample, collected with instantaneous luminosities up to 
~ 1.3 x 10 33 cm _2 s _1 . The trigger is fully efficient for the 4e, 4/x and 2e2/i channels. 

Monte Carlo (MC) samples for the SM Higgs boson signal and a variety of electroweak and QCD-induced 
SM background processes have been used for the optimization of the event selection prior to the analysis. 
They are also used for comparisons to the measurements, and the evaluation of acceptance corrections and 
systematic uncertainties. All production cross sections are re- weighted at least to next-to-leading-order (NLO), 
or beyond where possible. The general multi-purpose MC event generator PYTHIA |16j is used in conjunction 
with other event generators, for the showering, hadronization, decays and to add the underlying pp event. 
The events are processed with a detailed simulation of the CMS detector based on GEANT4 [TT]. The Higgs 
boson samples are generated with POWHEG [9 which incorporates NLO gluon fusion and weak-boson fusion. 
The events are re- weighted according to the most recent calculations [13] of the total cross section a(pp —> H) 
which comprises gluon and weak-boson fusion contributions. The total cross section is scaled by the £?(H — ► 4£) 
[5] [HI 03]. Interference in the 4e or 4/i channels is taken into account using the Prophecy4f [T3] generator tool 
for precision calculations. Di-boson production (WW, WZ, ZZ, Z7) is generated at leading order (LO) with 
PYTHIA. Interference effects involving the Higgs boson are neglected. An uncertainty of ±30 % is attributed the 
gluon-induced contribution. An inclusive Z+jets sample has been generated with MadGraph [5 , considering all 
(five) quark flavours for the initial state parton density functions (PDFs), and gluons, light (d,u, s) and heavy 
(c, b) quarks for the jets in the final state. A K-factor is used to correct to the NLO cross section. A top pair 
production sample is generated at NLO with POWHEG. 
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II. EVENT SELECTION 

The selection acts on loosely isolated (see below) lepton candidates, i.e. electrons within the geometrical 
acceptance of |?7 e | < 2.5 and p? > 7 GcV/c and muons satisfying |?7 M | < 2.4 and > 5 GeV/c. We require: 

1. First Z: a pair of lepton candidates of opposite charge and matching flavour (e + e~, /i + /i~) satisfying 

mi.2 > 60 GeV/c 2 , Pt,i > 20 GeV/c and p^,2 > 10 GeV/c; the pair with reconstructed mass closest to the 
nominal Z boson mass is retained and denoted 

2. Choice of the "best 4£": retain a second lepton pair, denoted Z 2 , among all the remaining £ + £~ com- 
binations with toz 2 > 12GeV/c 2 and such that the reconstructed four-lepton invariant mass satisfies 
rau > 100 GeV/c 2 . For the 4e and 4/z final states, at least three of the four combinations of opposite sign 
pairs must satisfy m« > 12GeV/c 2 . If more than one Z 2 combination satisfies all the criteria, the one 
built from leptons of highest px is chosen. 

3. Relative isolation for selected leptons: for any combination of two leptons i and j, irrespective of flavour 
or charge, the sum of the combined relative isolation Ri SO j + Ri SO ,i < 0.35. 

4. Impact parameter for selected leptons: the significance of the impact parameter to the event vertex, 
SIP 3 d, is required to satisfy |SIP 3 d = ^| < 4 for each lepton, where IP is the lepton impact parameter 
in three dimensions at the point of closest approach with respect to the primary interaction vertex, and 
(Tip the associated uncertainty. 

5. ZZ^*^ kinem atics: 60 < m Zl < 120 GeV/c and m% m < m Z2 < 120 GeV/c , where 

m min j g defied below. 

Requiring all selection criteria with m£"™ = 20 GeV/c 2 defines the baseline selection which is used in 
subsequent analysis, independent of the m H hypothesis. Requiring all selection criteria with m™ m = 60 GeV/c 2 
defines the high-mass selection which is used for the measurement of the ZZ cross section and for the Higgs 
boson searches at high mass (mn > 2 x m-£). 

The efficiencies for reconstruction, identification and selection for electrons and muons is measured with data 
by using a "tag-and-probe" technique based on the sample of inclusive single Z production events. The measure- 
ments have been performed in several ranges in \rj\ matching domains with uniform detector and reconstruction 
performances for the tracker and calorimeters, and in pr ranges from 7GcV/c to lOOGeV/c. 

The reconstructed leptons must also satisfy a very loose relative "track-only" isolation -Rfgo Ck < 0.7 where 
-Riso Ck = (1/VjO x ^so° k and Aj™ ck = J2iPr track- The sum runs over the tracks i within a cone of radius 

AR = VAT? 2 + A(j) 2 < 0.3 with axis along the lepton candidate direction, excluding a central inner "veto" 
region. For the leptons chosen to form the A£ system in event candidates, the combination of the tracker, 
ECAL and HCAL information is then used for isolation. The combined lepton R iso is calculated as Ri SO = 
(l/p l T ) x (A\™ ck + Af£ AL + Ago AL )- Thc ^f S o AL and ^So AL are sums ovcr thc e t from energy deposits in 
cells of the ECAL and HCAL respectively, and with geometrical centroids situated within a cone of radius 
AR — \J Aw 2 + A(j) 2 < 0.3. The measurements are performed for different \r]\ ranges covering the geometrical 
acceptance and px ranges up to 100GeV/ c. The combined isolation efficiencies measured with data using tag- 
and-probe is found to be above 99% everywhere for muons and between 94% and 99% for electrons. 

The overall signal detection efficiency for a A£ system within the geometrical acceptance is evaluated from MC 
to be rising from about 42% / 72% / 54% at m H = 190GcV/c 2 to about 59% / 82% / 71% at m H = 400GeV/c 2 
for the 4e / 4/i / 2e2 fi channels. A fit of the signal mass distribution as obtained from the MC shows a resolution 
for a Higgs boson mass hypothesis of 150 GcV/c in the 4e (4/x, 2e2^) of 2.7(1.6, 2.1) ± 0.1 GeV/c . Just above 
the 2 x m z threshold, for m H = 190GcV/c 2 , the resolution for 4e (4^, 2e2/i) is 3.5(2.5, 2.8) ± 0.1 GeV/c 2 . 

III. MEASUREMENT AND CONTROL OF THE BACKGROUND 

The small number of events observed precludes a precise evaluation of the background using mass sidebands. 
We therefore rely on other methods that use data to estimate the backgrounds and the associated systematic 
uncertainties. We choose a wide background control region outside the signal phase space which is populated 
by relaxing some selection criteria. 
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A. ZZ W Continuum 



Two different methods have been used to determine N^ pect for the ZZ^*) di-boson continuum, a normalization 
to the measured Z rate and an estimate from MC simulation. The former is our reference method for the present 
low integrated luminosities. A direct measurement from sidebands, e.g. excluding the signal region around a 
given Higgs boson mass hypothesis, is affected by large statistical uncertainties and will become useful only at 
higher integrated luminosities. 

The method of normalization to the measured Z rate relies on the measurement of inclusive single Z production 
which is used to predict the total ZZ rate within the acceptance defined by this analysis, making use of the ratio 
of the theoretical cross sections for Z and ZZ production, and of the ratio of the reconstruction and selection 
efficiencies for the 2£ and 4£ final states. The ZZ^*' — > 4£ control region is defined here by the observed single 
Z inclusive rate for Z — >■ ££. The ratio of acceptance between the control and the signal region is then given by 
combining the theoretical cross sections and the selection efficiencies as obtained from MC simulation: 

qq-tZZ-tM , ag^ZZ^U JLZ^-U 

ncr ne _ °NLO J g LO v fc MC /-i \ 

-"theory X %C ~ PP ^Z^2l X e Z^2i ' \ L > 

"nNLO £ MC 

The cross section for ZZ^ production at NLO through qq annihilation and gg fusion are calculated with 
MCFM, while the cross section for Z at NNLO is calculated with FEWZ 2.0. The theoretical uncertainties on 
the ratio are computed varying for each A£ and 2£ final states both the QCD renormalization and factorization 
scales (toz, tuz/2, 2toz). 

The efficiency cmc is defined as the ratio between events passing all the selection criteria and events with the 
cut at generator level (mu > 12 GeV/c 2 ). At the reconstruction level, the ZZ(*) events must fulfill all the event 
selection criteria. The selection of Z events follows from the first step of the selection. 

In a complementary method the estimation of the number of ZZ events in any given mass range [mi, 1112] is 
obtained directly from the absolute rate predicted by the MC model simulation. 

The number of events and relative uncertainties from ZZ^ 4£ predicted by normalization to the measured 
Z rate and by the MC model simulation for an integrated luminosity of 1.13 fb _1 in the signal region in a mass 
range from 100 to 600 GeV/c 2 with the baseline and high-mass selections are given in Table [ij 



TABLE I: Number of ZZ background events and relative uncertainties in the signal region in a mass range from 100 to 
600 GeV/c 2 , estimated from normalization to the measured Z rate and from Monte Carlo simulation. 





channel 


Normalization to Z rate 


MC model simulation 


baseline 


^yZZ->4e 


2.76 ± 0.18 


2.77 ± 0.26 






4.10 ± 0.27 


4.24 ± 0.39 




jyZZ->2e2fi 


6.72 ± 0.45 


6.85 ± 0.63 


high- mass 


^yZZ->4e 


2.50 ± 0.17 


2.52 ± 0.23 




jyZZ^4n 


3.55 ± 0.23 


3.66 ± 0.33 




jyZZ->2e2fj 


6.10 ± 0.40 


6.22 ± 0.58 



B. Reducible and Instrumental Backgrounds 

The Zbb/cc and tt control region is defined taking a set of events with a pair of identified leptons with opposite 
charge and matching flavour (e + e~, /.t + /i~) that form a good Zx, as requested for the signal selection. For the 
other pair of leptons the flavour, charge, and isolation requirements are removed. Finally the SIP3D impact 
parameter cuts are reversed on the two leptons by requiring ISIP3D | > 5. Inverting the SIP3D selection cuts for 
two leptons ensures a negligible Z+jcts contribution in the four-lepton background control region. 

To extract the number of tt and Zbb/cc events in the four-lepton signal region, we exploit the knowledge and 
distinct features of the SIP3D distribution. The SIP3D distributions for the Z2 leptons of the tt and Zbb/cc 
backgrounds are uniform and of similar shapes. This is in sharp contrast with the expected SIP3D distribution 
for the signal which is concentrated at low SIP3D and steeply falling with increasing SIP3D- 

In order to control the Z+jets background, a four-lepton background control region is obtained by relaxing 
the cuts on isolation and identification requirements for two additional leptons. 
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Normalized to the integrated luminosity, the number of events from tt, Zbb/cc and Z+jets expected in the 
signal region in a mass range from mi = 100 GeV/c 2 to to 2 = 600 GeV/c 2 is given in Table |h| 

TABLE II: Number of background events for baseline and high-mass event selections in the signal region in a m^l range 
from 100 to 600 GeV/c 2 , estimated from data. 





baseline 


high-mass 


j\jZbb/cc,tt^4t± 
jyZibb/ cc,tt— f 2e2/A 


0.01 ± 0.02 
0.01 ± 0.01 
0.02 ± 0.02 




jyZ+jcts-^4e 
jyZ+jets — >-4^i 
Z+jets— f2e2^i 


0.37 ± 0.07 
0.06 ± 0.01 
0.39 ± 0.07 


0.14± 0.06 
0.004 ± 0.004 
0.15 ± 0.06 



IV. SYSTEMATIC UNCERTAINTIES 



The main sources of systematic uncertainties on the expected yields are summarized in Table |III| Systematic 
uncertainties have been evaluated from data for the trigger efficiency as well as for effects from individual lep- 
ton reconstruction, identification and isolation efficiencies, and from energy-momentum calibration. Additional 
systematics come from the limited statistics in the background control regions which propagate to the back- 
ground evaluation in the signal region. All major background sources are derived from control regions, and the 
comparison of the data with the background expectation in the signal region is independent of the uncertainty 
on the LHC integrated luminosity for the data sample. This uncertainty enters in the calculation of a cross 
section limit through the normalization of the signal. Also given for completeness in Table |III| is an evaluation 
of the systematics uncertainties on the Higgs boson cross section and branching ratios. 

TABLE III: Summary of the magnitude of systematic uncertainties in percent. 



Luminosity 


6 


Trigger efficiency 


1.5 


Higgs cross section 


17-20 


Higgs B.R. 


2 


Lepton reco/ID eff. 


2-3 


Lepton isolation eff. 


2 


Electron energy scale 


3 



V. RESULTS 



The reconstructed four-lepton invariant mass distributions obtained in the 4e, 4/Lt, and 2e2/i channels with 
the baseline selection are shown in Fig. [I] and compared to expectations from the SM backgrounds. Three (six, 
six) event candidates are observed in 4e (4/i, 2e2^t) final states, satisfying the baseline selection. 

The reducible and instrumental backgrounds are very small or negligible. The number of events observed, 
as well as the background rates in the signal region wit hin a mass range from toi — 100 GeV/c 2 to rn-2 = 



600 GeV/c 2 , are reported for each final state in Table IV for the baseline and high-mass selections. 



The measured distribution is seen to be compatible with the expectation from SM continuum production of 
ZZ^*' pairs. We observe JV^ s scllnc = 15 events for the baseline selection, in good agreement with the expectation 
of 14.4 ± 0.6 events from SM background evaluation. Six of the events are below the kinematic threshold of 
two on-shell Z's (mjj < 180GeV/c 2 ), while 1.9 ± 0.1 background events are expected. The probability that the 
background fluctuates to the observed number of events is 1.3%. The events are not clustered, excluding the 
interpretation as the standard model Higgs boson. However we note that the six events form three mass pairs: 
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FIG. 1: Distribution of the four-lepton reconstructed mass for the baseline selection in the (a) 4e, (b) 4p, (c) 2e2/i, 
and (d) the sum of the 4£ channels. Points represent the data, shaded histograms represent the signal and background 
expectations. The samples correspond to an integrated luminosity of £. = 1.13 fb~ . 



two are close to each of the following masses 122, 142 and 165GeV/c 2 , respectively. We observe 7V^« hmass = 
8 events for the high-mass selection compared to an expectation of 12.4 ± 0.5 events from SM background 
evaluations. This corresponds to a 13% probability for the background to fluctuate to a value < j\f^ hniass _ 

No attempt to recover electromagnetic final state radiation prior to imposing the analysis selection criteria 
is made. However, a posteriori all selected events were inspected for the presence of direct photons that might 
originate by inner bremsstrahlung from leptons in the final state. 

The high-mass event selection and analysis which imposes the presence of two lepton pairs with invariant 
masses in the range 60 < mg+ e - < 120GeV/c 2 is used to provide a measurement of the total cross section 
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TABLE IV: Number of events observed, background and signal rates for each final state in a mass range from m\ — 
100 GeV/c 2 to m,2 = 600 GeV/c 2 both for the baseline and high-mass selections. For ZZ, Z+jets, ti and Zbb/cc the data 
driven estimations are used, for WZ the Monte Carlo estimation is used. 





4e 


baseline 
4 M 


2e2/j 


ZZ 
Z+jet 
Zbb/cc, ti 
WZ 


2.76 ± 0.18 
0.37 ± 0.07 
0.01 ± 0.02 
0.006 ± 0.006 


4.10 ± 0.27 
0.06 ± 0.01 
0.01 ± 0.01 
0.006 ± 0.006 


6.72 ± 0.45 
0.39 ± 0.07 
0.02 ± 0.02 
0.024 ± 0.012 


All background 


3.13 ±0.19 


4.17 ±0.27 


7.14 ± 0.46 


mji = 150 GeV/c 2 
m H = 200 GeV/c 2 
m H = 250 GeV/c 2 


0.368 ± 0.007 
0.816 ±0.011 
0.656 ± 0.009 


0.637 ±0.009 
1.161 ±0.014 
0.885 ±0.010 


0.996 ±0.011 
1.907 ±0.018 
1.533 ±0.013 


Observed 


3 


6 


6 






4e 


high-mass 
4 M 


2e2^ 


ZZ 
Z+jet 
Zbb/cc, ti 
WZ 


2.50 ± 0.17 
0.14± 0.06 


3.55 ± 0.23 
0.004 ± 0.004 


6.10 ± 0.40 
0.15 ± 0.06 

0.006 ±0.006 


All background 


2.64 ±0.18 


3.55 ±0.23 


6.25 ±0.40 


m H = 150 GeV/c 2 
m H = 200 GeV/c 2 
ma = 250 GeV/c 2 


0.018 ±0.002 
0.765 ±0.011 
0.621 ±0.009 


0.030 ±0.002 
1.080 ±0.013 
0.833 ±0.010 


0.045 ± 0.002 
1.801 ±0.017 
1.442 ±0.013 


Observed 
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a(pp — > TIL ± X) x /3(ZZ — > A£). This measured cross section is obtained via: 

0-{pp -> ZZ + X) X B(ZL -> 41) = ^ N ° hs ^ ~ ^backfrch) (2) 

A« x ezz^a x £ 

Where i c h means 4e, 2e2/j,; A is the acceptance of the detector, the efficiency ezz^u is the ratio between the 
ZZ events that survive the selection cuts and the generated ZZ events within the acceptance of the detector. 
The fiducial and kinematic acceptance is defined by the fraction of events with four final state leptons satisfying 
p T > 7(e), 5(/i) within |r7 e | < 2.5 and |t/"| < 2.4. The generated ZZ events arc the fraction of the PYTHIA MC 
sample used in the analysis with 60 < mz < 120GeV/c 2 . 

The total cross section for a pair of Z bosons in the mass range 60 < mz < 120 GeV/c 2 is found to be 

a(pp -> ZZ ± X) x S(ZZ -> At) = 17.5t^(stat.) ± 0.7(syst.) ± l.l(lumi.) fb . 

The measured cross section agrees within about 1.5 standard deviations with the expectation from the SM [7]. 

The exclusion limits for a SM-like Higgs Boson are computed for a large number of mass points in the mass 
range 110-600 GeV/c , and using the predicted signal and background shapes. The choice of the spacing between 
Higgs mass hypotheses is driven by cither detector resolution, or the natural width of the resonance, depending 
on which dominates. The signal shape is determined using 17 simulated samples covering the full mass range. 
The shapes for each simulated sample are fit using a function (/ ) obtained as a convolution between a Breit- 
Wigner-like probability density function to describe the theoretical resonance line shape, and a Crystal-Ball 
function to describe the detector effects. The parameters of the Crystal-Ball function are interpolated for the 
Higgs boson mass points where there is no simulated sample available. The mass shapes for the backgrounds 
are determined by fits to the MC samples, with full detector simulation, while the normalization is taken from 
the overall event yield estimates as described in previous sections. 

As a cross check the cut-and-count method has been used. To calculate the exclusion limits, the shapes are 
integrated in a given mass window to extract the signal and background yields. The mass window is chosen 
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FIG. 2: The mean expected and the observed upper limits at 95% C.L. on <r(pp — > H + X) x B(ZZ 4£) for a Higgs 
boson in the mass range 120-600 GeV/c 2 , for an integrated luminosity of 1.13 fb _1 using the CL S approach. The expected 
ratios for (a) the SM and for (b) the SM with a fourth-fermion family (SM4) are presented. The results are obtained 
using a shape analysis method. 



for each mass point to find the optimal expected limit, computed with a Bayesian approach with fiat prior. 
The resulting mass windows are asymmetric and their width accounts for the mass migrations from generated 
to reconstructed quantities, the mass resolution and the intrinsic width for a given central mass hypothesis. 
Similar exclusion limits results are obtained for the cut-and-count method. 

The observed and mean expected 95% CL upper limits on Higgs a(pp —> H + X) x B(ZZ — ► 4£) in a shape 
analysis method, obtained for Higgs masses in the range 110-600 GeV/c 2 are shown in Fig. [i] The limits are 
made using a CL S approach, for the expected ratios to the SM, and in the context of a SM extension by a 
sequential fourth family of fermions with very high masses (SM4) [T2 HH HH]. The bands represent the ler 
and 2cr probability intervals around the expected limit. The expected background yield is small hence the la 
range of expected outcomes includes pseudo-experiments with zero observed events. The lower edge of the 1<7 
band therefore corresponds already to the most stringent limit on the signal cross section, as fluctuations below 
that value are not possible. We account for systematic uncertainties in the form of nuisance parameters with 
a log-normal probability density function. The exclusion limits extend at high mass beyond the sensitivity of 
previous collider experiments. The expected limits reflect the dependence of the branching ratio £>(H — s- ZZ) 
on the Higgs boson mass. The worsening limits at high masses arise from the decreasing signal cross section. 
The reduced sensitivity around mn = 160GeV/c 2 and for the low masses arise from the very small H — > ZZ 
branching ratio in these regions. 

In the current study the exclusion limit is compared to the cross section for on-shell Higgs production and 
decay in the zero-width approximation, and acceptance estimates are obtained with Monte Carlo simulations 
that are based on ad-hoc Breit-Wigner distributions for describing the Higgs-boson propagation. 



VI. CONCLUSIONS 



The results of a search for a standard model Higgs boson produced in pp collisions at t/s — 7 TeV and 
decaying in ZZ'*' have been presented for the first time by the CMS Collaboration in the leptonic Z decay 
channel ZZ*-*-* — > i£, with I = e, fi. Simple sequential sets of lepton reconstruction, identification and isolation 
cuts and a set of kinematic cuts have been introduced to define a common baseline for the search at any Higgs 
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boson mass toh in the range 100 < mjj < 600GeV/c 2 . The instrumental background from Z+jets and the 
reducible backgrounds from Zbb and tt, with misidentified primary leptons, are shown to be negligible over 
most of the mass, with a small contamination remaining at low mass. 

Fifteen events are observed in the 2e2/x, 4e and 4^ channels for an integrated luminosity of 1.13 ± 0.07 fb -1 , 
while 14.4 ± 0.6 events are expected from standard model background processes. The distribution of events 
is compatible with the expectation from the standard model continuum production of Z boson pairs from qq 
annihilation and gg fusion. No clustering of events is observed in the measured mu mass spectrum. Six of the 
events are below the kinematic threshold of two on-shcll Z's (toh < 180GeV/c 2 ), while 1.9 ± 0.1 background 
events are expected. The probability that the background fluctuates to the observed number if events is 1.3%. 
Using the high-mass selection which contains eight events, a total cross section for a pair of Z bosons in the mass 
range 60 < to z < 120GeV/c 2 has been measured to be in agreement with the predicted value. Upper limits 
obtained at 95% CL on the cross section x branching ratio for a Higgs boson with standard model-like decays 
exclude cross sections from about one to two times the expected standard model cross section for masses in the 
range 180 < tor < 420GeV/c 2 . Upper limits obtained in the context the standard model with a fourth fermion 
family, exclude a Higgs boson with a mass in the ranges 138-162 GeV/c or 178-502 GeV/c at 95% CL. 
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